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Abstract 

The final state interactions of pions in decays K — ► ir^n ^ are 
considered using the methods of quantum mechanics. We show how 
to incorporate the electromagnetic effects in the amplitudes of these 
decays and to work out the relevant expressions valid above and below 
the two charged pions production threshold M c = 2m. The electro- 
magnetic corrections are given as evaluated in a potential model. 

During the last years essential progress has been achieved in tttc scattering 
lengths determination from experimental data [HE], [3]. The precise knowledge 
of these quantities is an important task, since at present the Chiral Pertur- 
bation Theory (ChPT) predicts their values with very high accuracy [I]. 
The high quality data on —>■ 7r ± 7r 7r° decays have been obtained recently 
in the NA48/2 experiment at CERN SPS [3J. Dependence of the decay rate 
on the invariant mass of neutral pions M 2 = (pi + P2) 2 reveals a prominent 
anomaly (cusp) at the threshold, relevant to the production of two charged 
pions Ml = 4m 2 . 

The theoretical investigations of the decays 

K± -> t^ttV, (1) 



K* -> t^tt+tT (2) 
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have been carried out for many years [SI El Ej.New experimental data of high 
accuracy [3] lead to great activity on this issue [HI El [TO] - 
As was explained by N. Cabibbo [8j,the cusp in the experimental decay dis- 
tribution is a result of the charge exchange scattering process 7r + 7r~ — > 7r°7r° 
in the decay (2). He proposed a simple re-scattering model,in which the am- 
plitude of decay (1) consists of two terms 

T = T + 2ika x T + (3) 

where T , T + are the "unperturbed" amplitudes for decays (1) and (2) re- 
spectively and k = ~V 'M 2 — 4m 2 is the momentum of the charged pion. 
The second term in (3) is proportional to the difference of scattering lengths 
o-x = (^o ~ fl2)/3 and flips from dispersive to absorptive at the threshold M c . 
As a result the decay probability under the threshold linearly depends on 
a x allowing this difference to be extracted from experimental data with high 
accuracy. 

The next important step was done in [9], where the amplitude T was ob- 
tained accounting the second order in scattering lengths terms using ana- 
lyticity and unitarity of the S matrix. The results of [9J were supported by 
rigorous approaches [TUHH], whose authors investigated the decay (1) in the 
framework of effective field theory and ChPT and calculate all re-scattering 
processes in two-loop approximation. The results from [pj were used in the 
fit of experimental data [3], allowing extraction of the difference ao — (i2 with 
high accuracy. 

Nevertheless the decay rate behavior near threshold cannot be provide solely 
by strong interaction of pions in the final state. As was widely discussed (see 
e.g. [121 EE]), the tiny discrepancy between theoretical predictions and ex- 
perimental data in the vicinity of threshold [3] is a result of disregarding 
electromagnetic effects in final state interaction. In view of importance of 
the knowledge of the scattering lengths with the most possible accuracy, the 
consideration of electromagnetic corrections in decay (1) becomes a pressing 
issue. 

Later on we discuss the problem of Coulomb interaction among charged pions 
using methods of non-relativistic quantum mechanics, which are completely 
suitable for the considered caseS- We obtain compact expressions for the am- 
plitude d of decay (1) with regard to the electromagnetic corrections, which 
are valid below and above the charged pions production threshold M c = 2m. 
Leaving the strict derivation for a separate publication, we shortly discuss 

lr The effects of radiation of real photons are beyond the scope of our consideration and 
will be treated elsewhere 

2 As the K~ decays are counterparts to the K + decays, they are not treated separately. 
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how to involve the electromagnetic effects in the considered problem and rel- 
evant modifications, which have to be done in the amplitude of decay (1). 
With the methods of non-relativistic quantum mechanics it can be shown 
that the result of N. Cabibbo [8] can be generalized accounting the nir scat- 
tering to all orders in scattering lengths 



T = T + 2ikf x T + , f x = a x /D, 

D = (l-ik 1 R 11 )(l-ik 2 R 22 ) + k 1 k 2 H 



(4) 



4m§; k 2 



k = \\fM 2 — 4m 2 are the momenta of neutral 



Here k\ - v , ,,. n . 

and charged pions respectively. The elements of the R matrix are real and 
can be expressed in isospin symmetry limit through the combinations of 
the scattering lengths [HI HD] a x = (cl 2 — a )/3;a o = (a + 2a 2 )/3;a± = 
(2a + a 2 )/6 corresponding to inelastic and elastic pion-pion scattering as 



R 



12 



V2a x ; Ru — a 00 ; R 



22 



2a, 



(5) 



The replacement a x f x has small numerical impact on the results of the 
previous calculations done according [HI [TO] in the dominant part of phase 
space, but, as we will see later, is very crucial for inclusion of the electro- 
magnetic interactions under the threshold,where formation of bound states 
(/T + 7r~ atoms) take place. 

The next step of our prescription is inclusion of electromagnetic effects in 
expression (4). The general receipt is known for many years (see for instance 
the textbook [H]) and implies replacement of charged pion momenta k by a 
logarithmic derivative of the pion wave function in the Coulomb potential at 
the boundary of the strong field r i.e. 



ik 



T 



dlog[G {kr) + iF {kr)\ 



dr 



(6) 



r=r 



Here Fq, Gq are the regular and irregular solutions of the Coulomb problem. 
In the region kr -C 1 where the electromagnetic effects can be significant 
the above replacement gives 



T 



ik — am 



log(-2i£r ) + 2C + i\) 1 



ima 
~2k 



Re(r) + ilm(r) 



Re{r) 
Im{r) 



-am 



\og(2kr ) + 2C + Re^[l- l -^\ 



,kA\ A = exp(^J|r(l + <)l, 



(7) 
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where C = 0.577, a = 1/137 are the Euler and fine structure constants, 
whereas ip is the digamma function. 

To go under the threshold, it is enough to make the common replacement 
k — > in in the above expression 



t = —k — am 



/ ma 

log(2«ro)+2C + V 1- — 



At K n = am I (2n), where n is an integer, r goes to infinity, which corresponds 
to Coulombic bound states in the considered approach. On the other hand, 
the product nf x defining the amplitude behavior under the threshold, remains 
finite due to dependence on r in the denominator D in expression (4). This 
explains why electromagnetic effects can be included only after summing up 
all terms of the infinite series in the perturbation expansion. 
The product nf x possesses a resonance structure placed at the positions 



AC? 



2m- -A 



am 



m 



— arctan A, A = am 



2(n-S)' 

«22 



1 + k\a\ x _ 



with the relevant width 



Airkia 



2 ^3 



m(l + klad) 



(9) 



(10) 



The physical reason of resonance origin is transparent. Due to the charge 
exchange process 7r + 7r~ — > 7r°7r° the Coulombic bound states of the 7r + 7r~ 
system (A 27T atoms) becomes unstable^. 

The considered effect of the creation of A 27T atoms in decay (1) is not the 
only contribution from electromagnetic interaction of pions. 
Outside the resonance region the Coulomb interaction leads to the essential 
difference between the r values calculated with electromagnetic corrections 
and without them. In particular, the nonzero contribution of the Coulomb 
corrections to the Re(r) above the threshold leads to the interference term 
in the decay rate provided by "direct" and "charge-exchange" contributions 
from (4). Thus above the threshold the interference is nonzero even at the 
lowest order in scattering lengths, unlike the original approach proposed by 
N. Cabibbo [8]. 

Further improvement of the theory consists in taking account of final state 



3 We do not discuss here the instability of excited states caused by their transition to 
the ground state, as this effect is very small and can safely be neglected. 
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interactions in the "direct" term from (4). This can be done by simple 
substitution 



T -> 7o(l + ifci/oo) 



(11) 





(12) 



(1 — z7cian)(l — TCI22) — ik\Ta\ 
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These higher order corrections to "direct" term are numerically small, but 
taking in account the precision of experimental data, have to be included in 
the fit procedure. 

To estimate the contribution from electromagnetic effects to the decay rate 



tude T is given by expression (4), while the amplitude T c taking account of 
electromagnetic effects is given by expressions (4,11,12) with relevant mod- 
ifications discussed above. The dotted line in Fig.l shows the contribution 
of electromagnetic effects without bound state corrections. The dashed line 
represents the same quantity, but with the corresponding average with the 
gaussian distribution [3] and the expected mass resolution near the threshold 
r.m.s. = 0.56 MeV. The solid line gives the contribution of all electromag- 
netic effects (bound states included) averaged as in the previous case with 
the gaussian distribution. From this plot one concludes that the essential 
contribution to the decay rate in the vicinity of the threshold comes from the 
electromagnetic interactions that do not lead to bound states. 
The developed approach allows one to take into account electromagnetic ef- 
fects in decay (1) and to estimate their impact on decay rate of the process 
under consideration. 

We are grateful to V. Kekelidze and J. Manjavidze who draw our attention 
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